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Abstract. We report a measurement of limb darkening of a solar-like star in the very high magnification microlensing event 
MOA 2002-BLG-33. A 15 hour deviation from the light curve profile expected for a single lens was monitored intensively in 
C/3 , V and I passbands by five telescopes spanning the globe. Our modelling of the light curve showed the lens to be a close binary 

system whose centre-of-mass passed almost directly in front of the source star. The source star was identified as an F8-G2 
main sequence turn-off star. The measured stellar profiles agree with current stellar atmosphere theory to within ~4% in two 
• »"H . passbands. The effective angular resolution of the measurements is <l/jas. These are the first limb darkening measurements 

' obtained by microlensing for a Solar-like star. 
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1. Introduction 

In gravitational microlensing events binary lenses produce 
bounded regions of high magnification, known as caustics, 
on the magnification map projected onto the source plane 
(Schneider & Weiss 1986 1. The steep magnification gradients 
associated with these caustic curves may be utilised to resolve 
the surfaces of background source stars as they move across 
them. Observations of this effect have been used to measure 
limb darkening in Galactic Bulge K giant stars (Albrow et al 
[T9991 Albrow et al l2000l Fields et al l2003l . a G/K sub-giant 
(Albrow et al 200 1 a I, and an A dwarf in the Small Magellanic 
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Cloud (reporting an angular resolution of a few nas, Afonso et 
al 2003. 



In this Letter we present observations of the very high mag- 
nification microlens event MOA 2002-BLG-33. In this event, 
the centre-of-mass of the close binary lensing system moved 
into near perfect alignment with a Solar-like source star. The 
resulting high magnification provided ideal conditions for mon- 
itoring the source as it transited the central caustic. Our obser- 
vations yielded the most precise limb darkening measurements 
obtained by microlensing for a non-giant star and the first such 
measurements for a Solar-like star other than the Sun itself. 



2 F. Abe et al.: Probing the atmosphere of a solar-like star by galactic microlensing at high magnification 



o 




HJD - 2450000 

Fig. 1. Observed light curve, best fitting model, and residuals from the model for the peak of MOA 2002-BLG-33. The total 
dataset consists of 473 measurements in the I band from Boy den 1.5 m telescope (South Africa), 186 points in Ve from EROS 
lm (La Silla), 260 points in I from MDM 2.4m (Kitt Peak), 480 points in I from MOA 0.6m (New Zealand), and 21 points in 
I and 10 points in V from Wise lm (Israel). The inset shows the separation of the light curves in the Ve and I bands caused by 
limb reddening. 



2. Light Curve of MOA 2002-BLG-33 

The event MOA 2002-BLG-33 was first alerted by the MOA 
collaboration on 2002 June 18 in its ongoing microlensing alert 
program (Bond et al 2001 1. About 48 hours before its time of 
peak amplification on 2002 July 4, MOA identified this event 
as a high magnification event in progress and a second alert 
was issued. Subsequently, an intensive campaign of follow-up 
observations was carried out using a number of telescopes cir- 
cling the globe (Fig. Photometry of the event was derived 
using difference imaging analysis applied to the CCD images. 
This procedure was carried out separately on each passband of 
each telescope resulting in six sets of time series photometry. 
These are shown in Fig. ^ where they have been normalised to 
the microlensing light curve using the best fitting model (de- 
scribed below). 

The M-shaped feature on the light curve results from the 
source star entering and exiting the central caustic curve pro- 
duced by a binary lens. The detailed shape of the light curve 
depends on the characteristics of the lens, the trajectory of the 



source star, and the radius and limb darkening profile of the 
source star. The parameters associated with the binary lens are 
the lens mass ratio M t /M2, the separation a of the two lenses, 
the Einstein radius crossing time f E , the source trajectory im- 
pact parameter, mo, with respect to the centre-of-mass of the 
lens, the time fo when the source is closest to the centre of mass 
of the lens, and the angle, <p, between the source trajectory and 
the lens components. 

The source star is parametrized by its radius r s and an ap- 
propriate limb darkening model. For this we used the square- 
root law (ClaretEOnn} 

j^- =1-q(1-cos0)-^(1- VEoi0) (1) 

where c A and d A are the wavelength dependent limb darkening 
coefficients. 

Model binary microlensing light curves were generated nu- 
merically using an inverse ray tracing technique, implemented 
on a cluster computer (Rattenbury et al. 2002 1. The inverse 
ray tracing technique naturally accommodates the finite size 
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Table 1. Microlensing parameters of MOA 2002-BLG-33. 
Here denotes the Einstein radius of the lens, as determined 
by its total mass, M\ + M2, Do L and Dos the distances to the 
lens and source stars. The errors are 2<x limits calculated as in 
Avni (tT976l . 



Parameter 




Value 


Lens mass ratio 


M1/M2 


0.54+0.20 


Lens separation 


a 


0.11±0.01 r E 


Einstein time 


fE 


50.7+1.0 day 


Position angle 


<P 


112°±2° 


Source radius 


r s 


(2.25±0.05)xl0- 3 r E D 0& /D 0L 


Impact parameter 


«0 


(4.29±0.09)xlO- 4 r E 


Centre time 


to 


HJD2452460.496+0.004 


Limb darkening 




0.41+0.04 




d> 


-0.08+0.65 




C V E 


0.49+0.02 




d' 


1.15+0.56 



Fig. 2. The passage of the source star of MOA 2002-BLG-33 
through the lens-caustic. The magnification map of the lens is 
represented in blue, with regions of higher magnification rep- 
resented by lighter shading. The caustic is the diamond shaped 
line of very high magnification. The red lines indicate the track 
of the source star relative to the caustic. The position angle of 
the lens binary system with respect to the source trajectory is 
denoted by (f>. The source star is shown to scale and in realistic 
colour. Its angular diameter would be 2-4yuas. 

of the source star and the limb darkening profiles across its 
face (Wambsganss I1997> . To find the best fitting model light 
curve, we implemented a x 1 minimisation technique using 
a Metropolis Hastings Markov chain Monte Carlo procedure 
(Ford 2003 1. This avoids the problem of traditional minimisa- 
tion techniques which can be fooled by false minima due to 
numerical noise in the simulated light curves. 

The best fitting theoretical curve is shown in Fig.^ A total 
of 1 1 physical parameters plus 6 flux normalisation parameters 
(one for each passband-telescope set) are required to describe 
the light curve 1 . The best fitting physical parameters are given 
in Table [0 The fitting procedure yielded correlated values for 
the limb darkening coefficients. In Tabled we report the coef- 
ficients (c\, d' A ), which are related to the coefficients in Eqn.^ 
through the phase space rotation c' A = ca cos 33°+dx sin 33° and 
d' A = d\ cos 33° - cx sin 33°. In this representation, the d' A term 
contributes a small correction to /^(^//^(O). The limb darken- 
ing parameters allow a colour map of the face of the source star 
to be drawn. This is depicted in Fig.[2]along with its passage 
through the caustic. 

Our best best fitting model light curve and the associated 
caustic correspond to a close binary lensing system. Dominik 

1 Blending parameters were not required because difference imag- 
ing analysis removes the effects of blending due to nearby neighbours 
and the lens itself. 



(1999) found that qualitatively similar caustic curves can be 
produced by lensing configurations that differ greatly. Caustics 
similar in shape to the one shown in Fig. [2] can be produced 
by lenses comprised of wide binary systems. Our best fitting 
wide binary solution for our observations of MOA 2003-BLG- 
33 yielded a mass ratio corresponding to that of a giant planet 
orbiting a stellar mass object. While this configuration qualita- 
tively produced a short duration M-shaped amplification pro- 
file, it could not reproduce the observed light curve profile at 
a level demanded by the precise, densely sampled photometry 
we obtained. Only the above close binary model gave an ac- 
ceptable fit to the data. 

Unfortunately, neither the rotation period of the binary lens, 
nor the plane of rotation, was able to be determined from the 
data. We estimated the contribution of binary rotation to our 
measurements of limb darkening as follows. For typical val- 
ues of the lens parameters, the rotation period lies in the range 
100-200 days. During the caustic crossing shown in Fig. [2] the 
lens is thus expected to rotate through 1 °-2°. Light curves were 
computed for rotating lenses with various periods of rotation, 
but with the plane of rotation held fixed and perpendicular to 
the lens-source axis, to maximise the effect. It was found that 
periods of rotation <150 days yielded light curves that were in- 
consistent with the data at a level greater than 2cr. Limb dark- 
ening coefficients, and stellar profiles, were then computed, in 
the static approximation, for a simulated dataset with a rotation 
period of 150 days. It was found that the computed stellar pro- 
files differed from the inputted profiles by approximately 2x 
the 2<x statistical uncertainty of the profiles computed without 
rotation. The effect of rotation was therefore allowed for, in the 
present work, by doubling the errors in the stellar profiles cal- 
culated in the static approximation. 

3. The surface of the source star 

We used colour measurements to identify the source star type 
in MOA 2003-BLG-33. The baseline (i.e. unmagnified) flux 
of the source star was determined by from the fluxes measured 
by the MDM and Wise telescopes at high magnification. The 
instrumental magnitudes were calibrated using the catalogue of 
Udalski et al ( 2002 1. The magnitudes thus determined in the V 
and I bands (Landolt system) were my = 19.94 + 0.05 and mi = 
17.96±0.07. Assuming a distance to the source star of 7±2 kpc, 
and absorption and reddening due to galactic dust with A\ = 
1.18 + 0.05 m ag and E(V-I) = 1.24 + 0.05 (Sumi etal. l20T)3l 
Udalski 2003 1, we derive its absolute magnitude M\ = 2.6 + 0.6 
and colour index V - I = 0.74 + 0.10. These values identify the 
source star as an F8-G2 main-sequence turn-off star with age 
approximately 3-10 Gy and temperature in the range 6200- 
5800K (Girardi et al. 2002 1, i.e. a solar-like star that has just 
commenced evolving away from the main sequence, or is just 
about to do so. 

Limb-darkening coefficients have been calculated for the 
ATLAS and PHOENIX models of stellar atmospheres for 
a wide variety of stars (Claret 12000b . The limb darken- 
ing coefficients, (c' A , d' A ), for the ATLAS model of a star 
with r e ff=6000K, surface gravity logg=4.0, and metallicity 
[M/H] = -0.3 are (0.363, 0.459) for the I-band and (0.461, 
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Fig. 3. Comparison of the measured profiles at visual (left) and 
infrared (right) wavelengths of the source star of MOA 2002- 
BLG-33, with predictions of the ATLAS and PHOENIX mod- 
els of stellar atmospheres (Claret 2000). 95% confidence limits 
are shown for the measured profiles. The x-axis variable, sin#, 
is the distance from the centre of the face of the star normalised 
to the star's radius. For the model atmospheres, the tempera- 
tures have been set at 5800K (lower curves) or 6200K (upper 
curves), the metallicity at [M/H] = -0.3 (Tiede et al. 1995 1, 
and the surface gravity log g = 4.0. 

0.394) for the V E -band 2 . These compare favourably with the 
measured values in in Table^that were derived from our obser- 
vations without including the effect of rotation described above. 
Fig- 13 shows the theoretical intensity profiles for the range of 
possible star types identified for MOA 2003-BLG-33, and the 
intensity profiles derived from the measurements with the ef- 
fect of rotation included. It is apparent that the ATLAS and 
PHOENIX profiles are confirmed at the 95% confidence level 
to within about 4%. The precision of the present measurements 
resulted from the fact that the magnification was high through- 
out the caustic crossing, and that it could therefore be mon- 
itored accurately everywhere, in particular at the end points 
which are most sensitive to limb darkening and reddening, as 
shown in Fig.^ The effective angular resolution of our obser- 
vations of MOA 2002-BLG-33 is < IRq/1 kpc, i.e. <lyuas. 

Our results appear to be in agreement with those expected 
by stellar atmosphere models for the source star type identi- 
fied here. Limb darkening measurements for another solar-like 
star, and consistent with those reported here, have been ob- 
tained by Deeg et al (2001 1 using a planetary transit technique. 
The only other observations by gravitational microlensing of 
limb darkening on a main sequence star were obtained for the 
A dwarf source star in MACHO 1997-SMC-l (Afonso et al 
2000 1. These were only partly in agreement with stellar atmo- 

2 For the V E -band, the average values given by Claret 1 2000 1 in the 
V and R bands were used, to allow for the extended red transmission 
of the EROS V filter, and also the pre-filtering effect of galactic dust. 



sphere theory (Claret 2000 1. Limb darkening measurements in 
agreement with theoretical expectations have been reported for 
a G/K subgiant source in OGLE 1999-BLG-23 (Albrow et al 
2001ai, and for the K giant source stars in MACHO 1997- 
BLG-28 (Albrow et al [T9^9l and in MACHO 1997-BLG- 
41 (Albrow et al 2000). On the other hand, recent results on 
source star in EROS 2000-BLG-5, also a K giant, appear in 
strong disagreement with stellar atmosphere theory (Afonso et 
al 2001 Fields et al 2003 1 whereas no such discrepancy was 
observed in H a equivalent width variations obtained in spec- 
troscopic observations of the same event (Albrow et al 2001b I. 

4. Conclusion and outlook 

We have shown how a high magnification microlensing event, 
with the lens being a close binary system, can be utilised to pro- 
file the atmospheres of a dwarf star. Such stars, even if nearby, 
cannot be resolved using conventional techniques. Thus the 
high magnification microlensing technique may prove useful 
until new, major programmes with the Very Large Telescope 
Interferometer and the Keck Interferometer become fully op- 
erational (Segransan et al. 2003 Domiciano de Souza et al. 
2003 1. We note that although some previous measurements us- 
ing microlensing techniques yielded possible evidence for sig- 
nificant departures from conventional stellar atmosphere the- 
ory, no such evidence was obtained in the present measure- 
ments. Indeed, as the star under study here was quite similar 
to the Sun, no such departure was expected. 
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